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Abstract 

This paper is concerned with a diffusion model of phase-field type, consisting 
of a parabolic system of two partial differential equations, interpreted as balances 
of microforces and microenergy, for two unknowns: the problem's order parameter 
p and the chemical potential each equation includes a viscosity term - respec- 
tively, e dtfM and 6 dtp - with e and 5 two positive parameters; the field equations are 
complemented by Neumann homogeneous boundary conditions and suitable initial 
conditions. In a recent paper [5], we proved that this problem is well-posed and in- 
vestigated the long-time behavior of its (e, 5)— solutions. Here we discuss the asymp- 
totic limit of the system as e tends to 0. We prove convergence of (e, (5)— solutions 
to the corresponding solutions for the case e = 0, whose long-time behavior we 
characterize; in the proofs, we employ compactness and monotonicity arguments. 
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1 Introduction 

The system we study was proposed for mathematical investigation in [S]; as to modehng 
issues, its most directly relevant antecedents are two papers by Fried & Gurtin [S] and 
Gurtin [9] , and a paper by one of us [11] . 

A nonstandard phase- field evolution problem. The initial/boundary- value prob- 
lem we dealt with in |5] consists in finding two phase fields, the chemical potential n and 
the order parameter p, such that 



e dtfi + 2pdtfi + fidtp — Afi = infix(0, +oo), (1.1) 

5dtp- Ap + f'{p) = fx in X (0,+cx)), (1.2) 

d^p = d^p = on r X (0, +oo), (1.3) 

p{-,0) = Po and p( ■ , 0) = po in fi, (1.4) 



where Q denotes a bounded domain of M.^ with (sufficiently) smooth boundary F, and /' 
stands for the derivative of a double- well potential /. This nonstandard phase-field model 
can be regarded as a variant of the classic Cahn-Hilliard system for diffusion-driven phase 
segregation by atom rearrangement: 

dtp~KAp = 0, p = -Ap + f'{p). (1.5) 

Note, in ( 11.11) . the unpleasant nonlinear terms involving time derivatives, and the fact that 
we have taken the mobility coefficient k > equal to 1. Moreover, recall that equations 
(II. 5p are customarily combined so as to obtain the well-known Cahn-Hilliard equation: 

dtp = KA{~Ap + f{p)). (1.6) 



Fried & Gurtin's generalization of Cahn-Hilliard equation. In [H [9] a broad 
generalization of (II. 6p was arrived at, with three measures: 



(i) by regarding the second of (11.51) as a balance of microforces: 

div^ + 7r + 7 = 0, (1.7) 

where the distance microforce per unit volume is split into an internal part vr and 
an external part 7, and the contact microforce per unit area of a surface oriented 
by its normal n is measured by | ■ n in terms of the microstress vector ^0 



ii) by interpreting the first equation of (11.51) as a balance law for the order parameter: 



dtp = — div h + a, (1. 
where the pair [h, a) is the inflow of p; 



^In [7], the microforce balance is stated under form of a principle of virtual powers for microscopic 
motions. 
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(iii) by requiring that the constitutive choices for tt, ^, h, and the free energy density ip, 
be consistent in the sense of Coleman and Noll [2] with a postulated "dissipation 
inequality that accomodates diffusion" : 

dt^J + (tt - fi)dtp - ^ ■ V{dtp) + h-Vfi<0 (1.9) 

(cf. equation (3.6) in [9]). 

In [9], the following set of constitutive prescriptions was shown to be consistent with (iii): 

7r(p,Vp,/i) = /i-9pZ?(p,Vp), > (1.10) 
|(p,Vp)=9vp^^(p,Vp) J 

together with 

h = -MVp, with M = M{p,Vp,p,Vfi), (1.11) 

provided the tensor-valued mobility mapping M satisfies the residual dissipation inequal- 
ity ^ 

Vp ■ M(p, Vp, p, Vp) Vp > 0. 

With the use of (11.71) . (11.81) . (ll.lOp . and (ll.liP j^. a general equation for diffusive phase 
segregation processes is arrived at: 

dtp = div (mv (dpipip, Vp) - div {d^p${p, Vp)) - 7)) + ^; 

in particular, the Cahn-Hilliard equation (II. 6p is obtained by taking the external distance 
microforce 7 and the order-parameter source term a identically null, and by choosing 

V^(p,Vp) = /(p) + ^|Vpp, M = Kl. (1.12) 

An alternative generalization of Cahn-Hilliard equation. In [11], a major modi- 
fication of Fried & Gurtin's approach to phase- segregation modeling was proposed. While 
the crucial step (i) was retained, both the order-parameter balance (II. 8p and the dissipa- 
tion inequality (II. 9p were dropped and replaced, respectively, by the microenergy balance 

dte = e + w, e:=— div/i + a, w := —n dtp + $ ■ V{dtp), (1-13) 

and the microentropy imbalance 

dtT] > — div h + a, h := p,h, a := p,a. (1-14) 

A further new feature was that the microentropy inflow {h,a) was deemed proportional 
to the microenergy inflow {h,a) through the chemical potential p, a positive field; con- 
sistently, the free energy was defined to be 



ip -.= 6- p, ^r], 



(1.15) 
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with the chemical potential playing the same role as the coldness in the deduction of the 
heat equationj^ 

Combination of fll.l3p -( [TTT5|) gives: 

dti^ < -Vdtif^-') +fx^'h-Vfx~7rdtp + $- V{dtp), (1.16) 

an inequality that replaces (11. 9p in restricting a la Coleman-Noll the possible constitutive 
choices. On taking all of the constitutive mappings delivering tt, ^, 77, and /i, dependent 
in principle on p, Vp, /i, Vp, and on choosing 

= ^^(p, Vp,p) = -pp + /(p) + ^|Vp|^ (1.17) 

compatibility with (I1.16P implies that we must have: 

f 7r(p, Vp, p) = -dp^ip, Vp, fi) = f^- f'{p), 1 

I ?(p, Vp, p) = 9vpV^(p, Vp, p) = Vp, I (1.18) 

I ^(p, Vp, /i) = /i^9/,^(p, Vp, /i)= -pV J 

together with 

/i(p, Vp, p, Vp) = H{p, Vp, p, Vp) Vp, Vp ■ -H"(p, Vp, p, Vp) Vp > 0. 

If we now choose for H the simplest expression H = k1 , implying a constant and isotropic 
mobility, and if we once again assume that the external distance microforce 7 and the 
source a are null, then, with the use of (11.181) and (11.151) . the microforce balance (I1.7P 
and the energy balance (11.130 become, respectively, 

Ap + p-/'(p) = (1.19) 

and 

2pdtp + pdtp- kAp = 0, (1.20) 
a nonlinear system for the unknowns p and p. 

Two-parameter regularization. Compare now the systems (I1.19l) - (ll.20p and (11.50 : 
(I1.19P and (ll.5P o are one and the same 'static' relations between p and p, whereas (ll.20p 
is rather different from (11.50 ]^. for more than one reason: 

(Rl) is linear, (OH]) is not; 

(R2) the time derivatives of p and p are both present in (ll.20p : 

(R3) in front of both dtp and dtp there are nonconstant factors that should remain non- 
negative during the evolution. 

^As much as absolute temperature is a macroscopic measure of microscopic agitation, its inverse - the 
coldness - measures microscopic quiet; likewise, as argued in the chemical potential can be seen as 
a macroscopic measure of microscopic organization. 
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Thus, the system fll.l9l) - fll.20p deserves a careful analysis. We must confess that we 
boldly attacked this problem as is, prompted to optimism by the successful outcome of 
a previous joint research effort [3l H] devoted to tackling the system of AUen-Cahn type 
one arrives at via the approach in [11] for no-diffusion phase-segregation processes. Un- 
fortunately, the evolution problem ruled by f ll.19p - p.20p turned out to be too difficult 
for us. Therefore, we decided to study its regularized version f ll.ip -f lT^ . which we ob- 
tained by introducing the extra terms e dt^ in f ll.20p and 6 dtp in ( ll.lQp . for small positive 
coefficients e and 6. Motivations for the introduction of such terms are proposed and 
discussed in [5]; interestingly, while the second can be interpreted as a dissipative part 
of the distance microforce, so far we have not been able to find a convincing physical 
interpretation for the first. But, our present study demonstrates - so we believe - that it 
can legitimately be regarded as an efficient mathematical device. 

Limit as the first parameter tends to 0. In [5], by assuming (as we did in [3lll]) 

that /' is the sum of a strictly increasing function /( with domain (0, 1) that is 
singular at the endpoints, and of a smooth bounded perturbation /2 (to allow for a 
double- or multi-well potential /), we first proved existence of a strong solution (/i,p) 
to fll.ip -( 1T^ satisfying fi > and < p < 1 almost everywhere in 1] x (0, +oo) (of 
course, we stipulated that the initial data meet same requirements in Q). Then, under 
some additional technical assumptions, we showed that the component is bounded, and 
so is /'(p); as a consequence, p stays away from the threshold values and 1. These 
boundedness properties are very useful in proving uniqueness of solutions. 

In some sense, passing to the limit as the regularizing parameters tend to zero is the 
challenging final aim of our research project. For the moment being, we are able to deal 
with e and to deduce, by a rather delicate asymptotic analysis, an existence theorem for 
the limit problem. Precisely, we let e tend to zero and show that any weak or weak star 
limit of any subsequence of solutions {fj,s,Pe) to f ll.ip -( !L4|) yields a solution (p,p) to the 
resulting limit problem, which is obtained by putting 5 = in f ll.ip - fll.4p and rewriting 
the corresponding first equation in the form 

2dt{fip) — Afi = ^dtp rather than as 2pdtp + pdtp — Ap = 0. (1-21) 

This we do because it is not clear to us from the structure of the system whether a suitably 
regular representation for dtp could be recovered in the limit, while we are able to show 
that the time derivative dt{pp) actually exists, at least in some generalized sensejfl 

Here, as we did in j5], we also deal with the long-time behavior of the system. We 
prove that each element (/i^, p^) of the cu-limit set for a certain trajectory is a steady state 
solution to f ll.ip - fll.4p : therefore, in particular, p^ is a constant (cf. f ll.2ip and f ll.3p )FI 

An outline of our paper is the following: in Section [21 we carefully state assumptions 
and results; Section [3] is devoted to the proof of the convergence theorem, so as to deduce 

■^In fact, (jl.ip has the equivalent formulation: 

dt{efi + 2fip) - = fidtp, 

which singles out the time derivative of the auxiliary variable e/i + 2/ip for e > 0. 

■^Note that the steady state problem associated with both cases £ > and e = is the same. 



6 



Asymptotic analysis for a nonstandard diffusion model 



the existence of solutions to the hmit problem; finally, in Section HI we develop our 
argument for the characterization of the w-limit. 



2 Assumptions and main results 

First of all, we assume f2 to be a bounded connected open set in with smooth bound- 
ary r and set, for convenience, 

V ■= H\n), H := L'^{n), and W := {v e H'^{n) : d^v = 0}. (2.1) 

We endow these spaces with their standard norms, for which we use self-explaining no- 
tation like II ■ II y. However, we write \\ -Wh for the norm in any power of H as well. The 
symbol ( ■ , ■ ) denotes the duality product between V*, the dual space of V, and V itself. 
Since Q is bounded and smooth, the embeddings W G V G H are compact. Moreover, 
since V is dense in H, we can identify H with a subspace of V* in the usual way, i.e., in 
order that {u, v) = {u, v)h, where ( ■ , ■)h denotes the inner product in H, holds for every 
u E H and v G V . Then, also the embedding H dV* is compact. 

As in [5], we assume that 

f = fi + /2, where /i, /2 : (0, 1) — M are functions satisfying: (2.2) 
fi is and convex, /2 is C^, and is bounded, (2.3) 
lim/((r) = — oo and lim/{(r) = -|-oo. (2.4) 

As to initial data, we start with the assumptions 

fiQ E V and /Uq > a.e. in f2; (2.5) 
PoEW, < Po < 1 in fi; and /'(po) e H. (2.6) 

The reader is referred to the forthcoming Remark 12.51 for weaker conditions. 

Since we aim to let e tend to zero, we stress the dependence of the solution found in [5] 
on the parameter e. In that paper, for any fixed T > 0, the following a priori regularity 
is required: 

jj, e H\0,T;H)r]L\0,T;W), (2.7) 
Pe e W^'°°{0, T; H) n i/^(0, T; V) H L°°(0, T; W), (2.8) 
/ie > a.e. in Qt, (2.9) 
< < 1 a.e. in Qt and /'(p^) G L°°(0, T; H), (2.10) 

where Qt := fi x (0, T). We note that homogeneous Neumann boundary conditions follow 
from fl2.7l) - fl2.8p . in view of the definition of W (see (12.11) ). Then, the e-problem is: 

[e + 2pe)dtPe + /^e dtPe - ^Pe = Or 

dt{eps + 2/iePe) - Ape = Pedtpe a.e. in Qt, (2-11) 

6dtpe - Ape + /'(pe) = p, a.e. in Qt, (2.12) 

Pe(0) = Po and Pe(0) = po a.e. in fi. (2-13) 
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We recall the existence result of [5]. 

Theorem 2.1. Let T G (0,+oo), and assume that (Q-dZlD and fl23|) - (l2:6|) are satis- 
fied. Then, there exists a pair {ji^, p^) satisfying (12. 7p -( 12TT0|) and solving problem f l2.1ip - 



As to uniqueness, the result in [5l Thm. 2.2] holds for solutions that, in addition to (12. 7p - 
(I2.10p . have certain properties that, in turn, are guaranteed whenever the initial data 
fulfil the following conditions, additional to (12. 5p and (12. 6p : 



(see Thm 2.3]). Within such a framework, since T > is arbitrary, the existence of 
a unique solution (/i£,pe) defined for every positive time was ensured, and its long-time 
behavior was studied. 

Here, our first concern is to construct a global solution {pi,,p^) to problem (12.111) - 
(I2.13P that satisfies (I2.7p - (l2.10p for every finite T, without assuming the just mentioned 
stronger conditions on the initial data. We cannot ensure uniqueness, of course. The 
corresponding result reads: 

Proposition 2.2. Assume that (I2.2p - (l2.4p and (12.51) - (12. 6p are fulfilled. Then, there exists 
a pair {^e, Ps) '■ [0, +oo) W x W satisfying (I2.7l) - (l2.10p and solving problem (12. lip - 
fl2T3|) for every T G (0, +oo). 

Starting from any family {{p^, Pe)}e>o of solutions of this type, we then let e tend to 
zero. To do this, we need to assume that 



in addition to (I2.5l) - (l2.6p . Under this assumption, we show that p^ is bounded away from 
zero and that (/i^, p^) tends to some pair (/i, p) as e \ in a suitable topology, at least for 
a subsequence. Moreover, we determine the limit problem solved by {p,p). The a priori 
regularity we require for [p, p) on every finite time interval is the following: 



(I2T3]). 



/io G L°°(f2), inf po > 0, and sup po < 1 



inf Po > 0, 



(2.14) 



pE L^iO,T;H)nL\0,T; V), 



(2 
(2 
(2 
(2 
(2 



15) 
16) 
17) 
18) 
19) 



p G H\0, T; H) n L°°(0, T; V) H ^^(0, T; W), 
p > a.e. in Qt and inf p > 0, 



p < 1 a.e. in Qt and /'(p) G ^^(0, T; H), 
pp eW^^P{0,T;V*) for some pG (l,+oo); 



the corresponding limit problem is: 




for every v E V and for a. a. t G (0, T), 
Sdtp - Ap + f'{p) = p a.e. in Qt, 
(pp)(0) = poPo and p(0) = po a.e. in Q. 



(2 
(2 
(2 



20) 
21) 
22) 
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Remark 2.3. The last integral in f l2.20p makes sense because V C L^{Q) and fidtp 
belongs at least to L\0,T- L^^^iQ)), as a consequence of fl2.15p -f r2.16p . Note that fl2.20p 
incorporates the homogeneous Neumann boundary condition for yU in a generalized sense. 
Moreover, note that (12.191) implies that fip is a continuous y*-valued function, so that 
the first equality in fl2.22p has a precise meaning. On the contrary, no continuity for fi is 
ensured at the moment. 

Here is our convergence result. 

Theorem 2.4. Assume that f EI^ -f l^ . fl23]) - fl^ . and (EH]) are satisfied, and let 
{(/^£) Pe)}ee(o,i) a family of solutions to problem fl2.lip - fl2.13p satisfying fl2.7p - fl2.10p . 
Then, there exists apair{n,p), satisfying fl2.15l) - fl2.19p and solving problem fl2.20p - fl2.22p 
for every T G (0,+oo), such that {fi^, Ps) converges to (/i,p) in a suitable topology)^ at 
least for a subsequence \ 0. 

Remark 2.5. The assumptions fl2.5p - fl2.6p are strong. In fact, while they are needed for 
Theorem 12.11 and Proposition 12. 2[ some of them will not play any role in the following, 
as it can be seen by looking at our a priori estimates. For instance, the last condition 
in (12. 6p will not be important, since just /(po) ^ L^{fl) will be used. Accordingly, one 
can prove a result similar to Theorem 12.41 but involving ^-approximations of less regular 
initial data that, this notwithstanding, satisfy fl2.5p - fl2.6p . Precisely, suppose we assume 
that 

HoeH, po>0 a.e. in fi, Po ^ V, f{po)eL^{Q), and infpo>0. (2.23) 

Then, it would be possible to construct ^-approximations (poe; Poe) of such initial data 
(a^O; Po) that satisfy fl2.5p - fl2.6p and whose norms of type (12.230 remain bounded as e \ : 
e.g., as to poe, one could take the solution to the elliptic equation 

g ^° ~ ^Pos + fiiPoe) = a.e. in fi, 
supplemented by homogeneous Neumann boundary conditions. 

Remark 2.6. Theorem 12.41 and the previous Remark offer us the possibility of defining 
and obtaining a weaker solution to problem fl2.lip - fl2.13p (that is, also for the case e > 0), 
if one writes equation (11.10 in the form (12.110 . To see that this solution is weaker than 
the one provided by Theorem 12.1] it suffices to compare (12. 150 - 02.160 with fl2.7p - fl2.8p . 
On the other hand, we can just assume fl2.23p and point out that in this approach one 
should consider fl219D . (KM . ((222]) with /xp replaced by (e/2)p + pp. 

Our final aim is to study the long-time behavior of any solution constructed according 
to Theorem 12.41 To this end, we introduce the cu-limit of the trajectory in a proper 
topology, and prove that every element of it coincides with a steady state. We set: 

w(/i,p) = {{p^,p^) E H XV : (p(tn),P(in)) (Pa;,Pa;) 

weakly in if x V for some sequence tn /" +C)o}. (2.24) 



^to be specified in the course of the proof given in the next section. 
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The above definition has a precise meaning, because the pointwise values of the {H x V)- 
valued function (/i, p) are well defined thanks to the continuity properties stated in our 
next result. 

Theorem 2.7. Assume that fra-fl^. fl^ - fl^ and fOij) are satisfied, and let 
(/i,p) : [0, +00) H X V be given by Theorem \ 2.4\ Then, (/i, p) is bounded, and its 



components /i and p are weakly and strongly continuous, respectively. In particular, the 
uj-limit (12.241) is nonempty. Moreover, every element of u{fi,p) coincides with a pair 
(/is, Ps) such that 

Ps is a nonnegative constant, 

ps^W, 0<p, <1, f'{ps)eH, and - Ap, + /'(p,) = p, a.e. mVL, (2.25) 
i.e., it coincides with a steady state. 



We stress that the above result does not necessarily hold for all possible solutions. 
Indeed, it only deals with solutions obtained as limits of solutions to the e-problem as 
e \ 0. We also observe that there is no reason for the function p^ mentioned in the 
statement to be a constant, since / is not required to be convex. 

The rest of the paper is organized as follows: in the next section, we prove both 
Proposition 12.21 and Theorem 12. 4t the proof of Theorem 12.71 will be given in the last 
section. 



3 Global solutions 

We first prove the existence of a global solution to the e-problem. The major part of the 
section is devoted to the proof of Theorem 12.41 and the subsequent existence of a global 
solution to the limit problem. 

Proof of Proposition [2T2l We imitate, with minor changes, the proof of Thm 2.1 in 
[S], where the final time T was fixed once and for all. Let e be fixed and, for notational 
conciseness, let the dependence on e be omitted. The main tool used in ^ was an approx- 
imation procedure using a time delay r = T/iV, for N a positive integer. Approximating 
T-problems were constructed and solved step by step in the time intervals /„ := [0,r;,r], 
n = 1, . . . , iV. It turned out that the resulting unique solution (p"^, p'^) coincided a poste- 
riori with the one obtained by glueing together solutions on the time steps [{n — l)r, nr], 
n = l,...,Ar. 

The necessary slight modification is the following. Take, e.g., t = 1/N, and solve the 
same problems as before step by step, now for every n > 1. This provides a global solution 
(p'^,p'^) to the approximating r-problem. Then, for every fixed T > 0, the argument of 
[5] applies, and a solution for (I2.1ip - (I2.13I) on [0, T] is constructed as the limit of the 
approximating solutions as r \ 0, at least for a subsequence. This holds, in particular, 
for T = 1, 2, . . . . Therefore, there is a subsequence ri „ \ such that the restriction of 
(p'^,p'^) to [0, 1] converges to a solution to problem (I2.11I) - (I2.13I) with T = 1. We denote 
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this solution by (/ii,pi). Now, take the restriction of (/i^,p^) to [0,2] with r = ri^„. 
Then, for the same reason, there is a subsequence {T2,n} of {'^i,n} such that the restriction 
we are considering converges to a solution {fi2, P2) to problem fl2.lip -f l2.13p with T = 2. 
However, as {t2,.„} is a subsequence of {ti „}, the restriction of {fJ^2, P2) to [0, 1] coincides 
with (/ii,pi). Proceeding inductively in this way, and then using a diagonal procedure, 
leads to a global solution to problem fl2.11l) - fl2.13|] . □ 

Preliminaries to the proof of Theorem 12.41 We begin by listing some of the tools 
we shall use. First of all, the well-known continuous embedding, with the related Sobolev 
inequality, holds in our 3-dimensional case: 

W^'P{n) C L'^in) and \\v\\L,{n) < Cp\\v\\w^.p(^n) for every v G W^'P{n), (3.1) 

provided that 1 < p < 3 and 1 < q < p* '■= , (3-2) 

3 — p 

with the constant Cp in (13. ip depending only on Q and p; moreover, 

the embedding W^'^in) C L\n) is compact if l<q<p*. (3.3) 

In particular, V C L'^{Q) for 1 < q < Q, and 

||'y||L'3(n) < C'll'^lly for every v and 1 < g < 6, (3.4) 

where C depends only on Q and the embedding V C L^(f2) is compact if g < 6. Further- 
more (see, e.g., P formula (3.2), p. 8]), we have the continuous embedding 

L°°(0, T; H) n L\0, T; V) C L^'^^^iQx) 

and the related inequality 

II^IIlio/3(Qt) — C'T||'y||l,°°(0,T;H)nL2(0,r;y) 

for every v E L~(0, T; H) D ^^(0, T; V^), (3.5) 

where Ct depends on Q and T. In our proof, we shall make use also of the well-known 
Holder inequality, mainly in the form 



n 



\vi - ■ ■Vn\\LP{o,T;Li{n)) < Yl\\vi\\ Ln (o,T■,L1^ (n)) for Vi E L^' (0, T; L"' (fi) ) , i = 1, 

i=l 

1 " 1 1 " 1 

provided that p, q,Pi, qi E [1, +oo], - = — and - = — . 



Remark 3.1. To avoid a cumbersome notation, the lowercase letter c stands for different 
constants, each of which may depend on one or another of the data involved in our current 
statement and on the coefficient 6, but never depends either on e or on the final time T; 
consequently, the relative estimates continue to hold when we discuss both the system's 
asymptotic limit as e tends to and its long-time behavior. Moreover, a notation like c^r 
signals that that constant has an additional dependence on the parameter a. Hence, 
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the meaning of c and may change from hne to hne, and even in the same chain of 
inequahties. On the contrary, we use the uppercase letter C for precise constants we are 
going to refer to after their introduction, such as Cp in (13. ip or Ct in (13 ■4p . Finally, in 
order to lighten our notation, we do not write the subscript e in performing our a priori 
estimates until each estimate is completely proved; the same we do for the auxiliary 
function 

Ue := epe + (3.6) 

Next, we prove that is bounded away from zero uniformly with respect to e. Such 
a result is essentially known from the proof of Thm 2.3], among other properties there 
established for a fixed e. Nevertheless, we prefer to repeat the proof here, in order to 
make sure that the constructed lower bound is in fact independent of e, and that just the 
additional assumption (I2.14p is used. 

Lemma 3.2. There exists some G (0, 1) such that > r^, a.e. for every e G (0, 1). 

Proof. We set for convenience := inf po and M := sup^gj-o^i) |/2('")| and recall that 
p* > by (I2.14p . Thus, owing to (12. 4p . we can choose r* G (0,p*] such that /{(r^,) < — M. 
Then, we test (I2.2ip by — (pe — r*)~ and integrate over Q x (0,t) where t G (0,T) is 
arbitrary. By omitting the subscript e for simplicity, we have 

^ / l(p-r*)-(t)r+ fl |V(p-n)"r- fl {f[{p)-f[{r.)){p-r^)- 
^ Jn Jo Jn Jo Jn 

= U lip-r.rm'- fl p(p-r.)-+ fl {n{r.)-m){p-n)-. 

^ Jn Jo Jn Jo Jn 

Every term on the left-hand side is nonnegative; in the right-hand side, the first term 
vanishes, because po > r^,, and the other two are nonpositive, because p > and /{(r*) — 
/2(p) ^ f[{r^:) + M < 0. Hence, (p — =0, and the assertion is proved. □ 

Proof of Theorem 12.41 Our proof will proceed as follows. For a fixed finite final 
time T, we shall perform a number of a priori estimates and use well-known compactness 
results to prove that, as e tends to 0, the solution (pe,pe) to the e-problem (I2.1ip -( l2.13p 
we are considering converges to a solution (p, p) to problem (I2.20p -( l2.22p . at least for a 
subsequence En \ 0; in particular, this holds for T = 1, 2, . . . . Having established this 
result, we shall be able to argue as in the proof of Proposition 12.21 Indeed, a diagonal 
procedure provides a subsequence e„ \ such that {pe,Pe) with e = En converges to a 
global solution (p, p) to problem (I2.20I) - (I2.22I) defined in the whole of [0, +oo). Therefore, 
just the case of a fixed final time T has to be considered. 

First a priori estimate. We test (12.111) (e.g., the equation within square brackets) 
by pe, and integrate over Q x (0,t), for an arbitrary t G (0,T). We obtain 

fl dJlp' + pp')+ fl |Vpp = 0, 
JoJn ^ JoJn 

whence 

^ / Ht)\'+ I {Pf^'m+ I I \Vfi\' = ^-WpoWh + Wpot^lhHn) < c. 
^ Jn Jn Jo Jn ^ 
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Since p/i^ > r^,/i^ thanks to Lemma 13. 2[ we immediately deduce that 

\\Pe\\L^{S),T;H) + || V/i^ || L2(o,T;i/) < C. (3.7) 

Second a priori estimate. We test fl2.12p by dtPe, and use the second of (12.111) in 
order to compute the right-hand side we get; we also recall (13.61) . For t G (0,T), we 
obtain: 

5 fl \dtp?+\ I ivp(t)r-^ / ivpor+ / /(p(t))- / /(po) 

Jo Jn Jn ^ Jn Jn Jn 

nl^dtp = j dt{ep + 2pp) - / / A/i 
_ _j Jo Jn Jo Jn 

ndtu = / u{t) - / {epo + 2poPo) < 3 / fi{t) + c. 
. _| Jn Jn Jn 

Since (13. 7p holds and / is bounded from below, we easily infer that 

\\dtPe\\L^O,T;H) + W"^ Pe\\ {0,T;H) + \\ f iPe)\\ {0,T;L^n)) < C. (3.8) 

Moreover, because < p^ < 1 a.e. in Qt for every 5 G (0, 1), we conclude that 

\\Pe\\H^O,T;H)nL°°{0,T;V) < ^T- (3.9) 

Third a priori estimate. Taking into account (I3.7l) -( !3l9l) and the boundedness of /g, 
we see that (I2.12p yields 

\\-^Pe + f[{pe)\\L^O,T;H) = Wf^e " dtps - /a (Pe) ||l2(0,T;H) < Ct • 

By a standard argument (test formally by f[{pi;), for instance) and elliptic regularity, we 
conclude that 

\\f[{Pe)\\L^{0,T;H) + || Pe || (0,T;W) < Ct ■ (3.10) 

First conclusions. The above estimates allow us to use standard weak and weak star 
compactness results. Thus, a triplet (p, p, (p) exists such that 

Pe^/i weakly star in L°°(0,T;i7) nL2(0,T;\/) (3.11) 
Pe^p weakly star in 1/1(0, T; if) nL°°(0,T;\/) nL2(0,T;iy) (3.12) 
f[{pe)-^ip weakly in L2(0,T;i7), (3.13) 

at least for a subsequence En \ 00 We note that p > and that p > r^, a.e. in Qt (the 
former inequality holds because Pe > for every e, the latter by Lemma [3. 2p . Moreover, 
by (I3.12P and the compact embedding V C U\VL) for p < 6, we infer that 

pe^ p Strongly in C°([0, T]; L^{VL)) for every p < 6, (3.14) 



^Incidentally, we anticipate that the convergence results stated below will hold only for suitable sub- 
sequences. Nevertheless, we will not mention such a detail. 
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thanks to [T^l Sect. 8, Cor. 4]. Hence, /2(Pe) converges to /2(p) in a suitable topology, 
for instance, strongly in L^(0,T;if), since is Lipschitz continuous. In particular, we 
deduce that 

5dtp - Ap + if + f^ip) = p a.e. in Qt- 

Furthermore, invoking both (13.140 and (I3.13p . and using a standard monotonicity tech- 
nique (see, e.g., [H Lemma 1.3, p. 42]), we conclude that 

< p < 1 and if = f[{p) a.e. in Qt- 

Finally, fl3.14p implies that Pe(0) converges to p(0) strongly in H, whence p(0) = po- 

In summary, so far we have proved fl2.15p - fl2.18p . (12.210 . and the second condition 
in fl2:22|l . It remains for us to show fl2:T9|) . fl2:20D . and the first condition in (K22^ . For 
this purpose, further arguments are needed. 

Fourth a priori estimate. We recall (13.60 . and we notice that the second (12. lip reads: 

dtUe = PedtPe + Ap^ . (3.15) 

Moreover, p^ satisfies homogeneous Neumann boundary conditions, since it is VT-valued. 
Therefore, we have 

dtUe v = [ [ PedtPe V- ! / Vp^ " for every v G ^^(0, T; V^). (3.16) 
Jn Jo Jn Jo Jn 

From (13.160 . we derive an estimate for S^Ug as a \^*-valued function, in the framework of 
the Hilbert triplet {V, H,V*). We treat the integrals on the right-hand side individually 
(for a while, we omit the subscript e in order to simplify the notation). 

Assume that v G L^(0,T; V). Then the Holder inequality and the Sobolev inequality 
fra with q = 5 yield: 

f I fldtpv <\\p\\L^o/^Q^)\\dtp\\L^Qr)\M\L^Qr) 

Jn 

< c||p||iio/:i(Qy) \\dtp\\L^{Qj,) \\v\\l5{0,T;V) ■ 

On the other hand, inequality (13. 5p holds. Therefore, on taking into account (13. 7p 
and fj3.90 . we conclude that 

T 

pdtpv < CT\\v\\Li(Q^T;V) for every v G L^(0,T; V). (3-17) 

Jn 

Next, we consider the second term on the right-hand side of (13.160 . By assuming again 
that V G L^(0,T; V), and invoking (13.70 once more, we immediately find that 





T 




Vp-Vv 



< 



\p\\l'^{0,T;V) \\v\\l2{0,T;V) < Ct\\v\\ l5 (0,T;V) ■ 



'0 Jn 

Combining this estimate with (13.170 and fl3.16p . we obtain that 

\{dtu,v)\ < cr||f ||L5(o,r;y) for every v G L^{0,T]V). 
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In other words, we have that 

\\dtUe\\L5/i{Q^T-y*) < ct- (3.18) 

Consequence. From (i) the strong convergence (13.141) with p = 4, (ii) the weak conver- 
gence /ij — !■ /i in L^(0, T; L^(r2)) imphed by (13. lip , and (iii) the Sobolev inequahty (13. 4 p 
with g = 4, we infer that 

/iePe — /xp Weakly in L^(0, T; if ), whence — )■ 2pp weakly in L^(0, T; ff), (3.19) 

since ep^ — )■ strongly in L^(0,T;y), by (13. 7p . Hence, accounting for (I3.18p . we con- 
clude that 

dtUe -> 29t(pp) weakly in L^'^{Q,T-V*), whence -> 2pp weakly in W^^^/^{Q,T-V*), 

(3.20) 

so that (I2.19P holds with p = 5/4. Moreover, (I3.20p also implies that converges to 
2pp weakly in C°([0,T]; V*); in particular, ^^(0) (2/ip)(0) weakly in V*. On the other 
hand, Ue(0) = epo + 2popo converges to 2popo, e.g., strongly in H. Thus, the Cauchy 
condition for pp in (I2.22p follows. 

In order to prove (I2.20p . one can try to let e \ in (I3.16P first, then to get rid of 
time integration. But, a trouble arises in dealing with the first term on the right-hand 
side, since, for the moment being, both p^ and dtPe are just weakly convergent. Hence, 
we have to prepare a new tool. 

Fifth a priori estimate. We want to find a bound for Vu^, i.e., for the partial deriva- 
tives -DjUe, z = 1, 2, 3. As usual, we omit the subscript e for a while. We have: 

\Diu\ = \eDifi + pA/i + ^J'Dip\ < 2\Difi\ + p| Ap|- 

Now, on taking (13. 7p into account, we see that Difi is bounded in L^(0,T; H), while p is 
bounded in L'^{0,T; L^{Q)) thanks to the Sobolev inequality (13. 4p . On the other hand, 
(13.90 provides a bound for Ap in L°°{0,T; H). Hence, using Holder inequality, we see 
that the product pAp is bounded in L^(0,T; L^/^(f2)). Therefore, we conclude that 

ll^e||L2(0,T;Vyi.3/2(Q)) < Ct ■ (3-21) 

Consequence. As (I3.19P holds, from (13.2 ip we infer that 

Ue 2pp weakly in L^{0,T; W^'^/^{n)). (3.22) 

Now, we observe that the embedding W^'^^'^{il) C L'^{il) is compact for every g < 3, 
by (13. 3p . On the other hand, (I3.20p holds. By using the Aubin-Lions lemma (see, e.g., 
[TUl Thm. 5.1, p. 58]), we deduce the strong convergence 

Us -> 2pp strongly in L'^{0, T; Li{Q)) for every g < 3. (3.23) 

We stress that, in particular, Ue —J- 2pp strongly in L^(0,T; H). 
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Lemma 3.3. The strong convergence fi^ ^ fi holds in L^(0,T; H). 



Proof. We set u := 2/ip and argue a.e. in Qt for a while. Thanks to Lemma [3. 2 [ we have 



U 



2pUs — eu — 2p^u 



2p{e + 2p, 



< 



e\u\ + 2\pUi; — Psu\ 



e + 2p, 2p 
On the other hand, we have 

\pUe - PeU\ < \p\ \Ue - u\ + \u\ \p - Pe\ < \Ue - u\ + 2p \p - p^\. 

By combining these inequahties, we deduce that 

WfJ-e - fJ'\\L^{0,T;H) 

< C (£:||M||i2(o,T;/f) + Ike - 'u||L2(0,T;/f) + ||/^||L2(0,T;L4(n))||P - Pe||cO([0,T];L4(Q))) • (3.24) 



The first two terms on the right-hand side tend to zero as e \ 0, by (13.231) : as to the last 
term, it suffices to recall f EJf^ . flXTI]) . and ([S3D for g = 4. □ 

Conclusion. The strong convergence guaranteed by Lemma together with the weak 
convergence dtPe — )■ dtp in L^(0,T; H) given by 03.121] . imply that 

Pe dtPe P dtp weakly in (Qt). 



On the other hand, fl3.20p and (13. lip hold. Hence, by letting e \ in (I3.16p . we easily 
obtain that 



{dtM{t),z{t))dt 



pdtp z — I V/i ■ Vz 




Jn 




(3.25) 



Jn 



for every z G L^O,T; V)nL'=^{QT). Now, take any v eVnL'^{n) and any C e L°°(0,T). 
Then the function z : t C{t)v is admissible in (13. 25 p . and a standard argument yields 



2{dt{pp)it),v) + / Vp{t)-Vv= / p{t)dtp{t)v, 
Jn Ja 

for every v ^ V H L°°{Q) and for a. a. t G (0,T). Now, we note that, for a.a. t G (0,T), 
each term in the above equation defines an element of V*. This is clear as far the left- 
hand side is concerned, since dt{pp) is l^*-valued and p is valued. For the remaining 
term, we recall Remark 12. 3[ On the other hand, V fl L°°{Q) is dense in V. Therefore, 
(I2.20p follows, and the proof is complete. □ 



4 Long-time behavior 

This section is devoted to proving Theorem 12. 7[ We first derive some a priori estimates, 
then we prove the continuity property announced in the statement; finally, we characterize 
the cj-limit. 
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A priori estimates. We recall Lemma [3.21 and the a priori estimates (13. 7p and (13. 8p . 
which involve constants that do not depend on the final time. Hence, we immediately 
obtain that 

p{t) > n and \\fi{t)\\H + \\^p{t)\\H < c for a.a. t > 0. (4.1) 

Recalling that < p < 1, we see, in particular, that (/i, p) is a bounded {H x V^)-valued 
function, as stated. For the same reason, we also deduce that 

\\M{t)\\H<c fora.a. t>0. (4.2) 
Moreover, the estimates (13.71) and (13. 8p also yield the bounds 

||Vyu||L2(o,r;/^) + \\dtp\\L2{o,T;H) < c for every T > 0, 
and we conclude that 

/•OO /» /"OO /» 

/ / |V/x|^ < +00 and / / l^tpl^ < +oo. (4.3) 
Jo Jn Jo Jn 

Strong and weak continuity. As far as p is concerned, we have p G H^{0,T; H) fl 
L^(0,T; W) for every T < +oo by (I3.12p . Since the embedding 

H\0, T; H) n L2(0, T; W) C C°([0, T]; V) 

holds, we immediately deduce that p is a strongly continuous V"-valued function. The 
weak continuity of /i is less obvious: we prove it by using the following well-known tool, 
whose proof is left as an exercise to the reader. 

Proposition 4.1. Let Z be a H aus dor ff topological space, and let Z he a reflexive Banach 
space such that Z <zZ, where the embedding is continuous with respect to the weak topology 
of Z . Assume that z : [0,T] ^ Z is continuous and that z{t) G Z and \\z{t)\\z < M for 
a.a. t G (0,T) for some constant M. Then z is Z-valued, i.e., z(t) belongs to Z for 
every t G [0,T], and is continuous with respect to the weak topology of Z . Moreover, 
\\z{t)\\z < M for every te [0,T]. 

In our case, we argue on any fixed finite time interval [0, T] and apply Proposition 14. II 
twice, first with Z = V*, with either the weak or the strong topology, then with Z = L^{Q), 
endowed with the weak topology. We set: 

u := 2/ip, (4.4) 

in order to agree with ([XnD, and we recall that u G W^'^/^{0,T;V*), by flX^ : m par- 
ticular, u G C°([0,T]; V"*). On the other hand, we have proved (14. 2p . We conclude that 
u(t) G H for every t G [0,T] and that u is continuous with respect to the weak topology 
of H. Besides, p is strongly continuous as an //-valued function, and the first condition 
in (14. ip holds. Hence, the same is true for 1/p, and we infer that p = u/{2p) is a weakly 
continuous L^(r2)-valued function. Now, we recall the estimate of p given by (14. ip and 
conclude that p is weakly continuous as an //-valued function. 
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Conclusion. It remains for us to sliow tliat every element of the w-limit is a steady 
state. To this end, we pick any {fiuj, Puj) G u{fi,p) and consider a corresponding sequence 
tn /• +CXD, as given by definition f l2.24p . We set: 

/i.„(t) := + t), pn(t) := p(tn + t) , and Un{t) := u{t + tn) , for t > 0, (4.5) 

and study the sequence {{pn, Pn)} on a fixed finite time interval [0, T] by using Un as well. 
From (14. ip . (14. 3p . and weak star compactness, we immediately deduce that 

P^n poo weakly star in L°°(0, T; H), Pn ^ Poo weakly star in L°°(0, T; V^), 
\Vpn\ and dtPn strongly in L^(0,T; H), 

at least for a subsequence. It follows that p^o is space- and poo time- independent. Thus, 
we can write pooif) = Ps for a.a. t G (0,T) for some pa G V. On the other hand, we can 
reproduce the estimates (13.101) . (13.181) . and (I3.2ip . on the time interval [t„,t„, -|-T] instead 
of [0,T]. We obtain: 

where ct does not depend on n. This means that 

l|/l(Pn)IU2(0,T;H) + ||Pri||L2(0,T;W) + l^n II Tyi.5/4(0,r;y*)nL2 (0,T;VFi'3/2(f7)) < ■ (4.6) 

Thus, Poo £ -^^(0,T;iy), i.e., G W. Moreover, to derive a strong convergence for p„, 
in C''([0, T]; iJ), we can argue as in the previous section. This allows us to ensure that 
/2(pn) converges to /2(Poo) strongly in L^(0,T;if) and that the weak limit of /((p„) in 
L^(0,T; if), given by weak compactness, is f[{poo)- All this yields that 

< Ps < 1 and - Aps + f'{ps) = Poo a.e. in Qt, 

and we deduce that poo is even time-independent. Thus, Poo{x, t) = ps for a.a. (x, t) G 
for some constant ps- Furthermore, ps is nonnegative, since Pn > for every n. This 
concludes the proof that (ps,Ps) is a steady state. 

Lastly, we show that {ps,Ps) coincides with (p^,p^). Because Pn P strongly in 
C°([0, T]; if), we see that Pn(0) converges to Poo(O) = Ps strongly in H. On the other 
hand, by assumption p„(0) = p(t„) converges to p^ weakly in V. Hence, Ps = Puj- 
The corresponding argument for p^ and bit more involved. We remind that the 

embedding W^'^^'^{il) C if is compact. Hence, from (14. 6 p and the Aubin-Lions lemma, 
we conclude that there is some Uoo such that 

Un Uoo weakly in W^'^/\0, T; V*) n ^^(0, T; W^^^^'^iQ)), (4.7) 
whence Un ^ u^o strongly in L^(0, T; if). (4.8) 

On the other hand, the strong convergence p„ — > poo in C°([0, T]; if ) and the uniform 
inequality pn > r^, imply the strong convergence l/p„ 1/poo in C''([0, T]; if ). We infer 
that 

= ^ ^ ^ Strongly in L\0, T; L\n)). 
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Since fin — >■ /ioo weakly star in L°°{0,T; H), we conclude that Moo/(2poo) = A^oo, i-e., that 
Uooit) = for a.a. t G (0,T). Next, the first weak convergence (14.71) also implies weak 

convergence in C^{[0,T];V*)] in particular, Un{0) converges to Uoo(O) = 2iiap^ weakly 
in V* . On the other hand, by assumption /i(t„) — ?■ /i^ weakly in H and, due to the 
already mentioned strong convergence p„ — )■ poo, p(^n) Pu strongly in H. We infer that 
UniS^) = '^f^{tn)p{'tn) couvergcs to 2ji^p^ weakly in L^iVt). By comparison, we conclude 
that 2psPuj = ^pi^pi^, i.e., that ps = Poj- CH 
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